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The optimized linear (^-expansion is applied to multi-field 0{N\) x 0{N2) 
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for parameters values for inverse symmetry breaking (or symmetry nonrestora- 
tion) at high temperatures. Our results are compared with the ones obtained 
by the one-loop perturbative approximation, the gap equation solutions and 
the renormalization group approach, showing good agreement with the lat- 
ter method. Apart from strongly supporting inverse symmetry breaking (or 
symmetry nonrestoration) , our results reveal the possibility of other high tem- 
perature symmetry breaking patterns for which the last term in the breaking 
sequence is 0{Ni - 1) x 0{N2 - 1). 
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I. INTRODUCTION 



The possibility that symmetries may be broken (or remain broken) at high temperatures 
is not new This phenomenon is usually called inverse symmetry breaking (ISB) (or 
symmetry nonrestoration (SNR)). The idea of ISB (or SNR) is per se a very interesting one 
due to its possible implementation in realistic particle physics models and its consequences 
in the context of high temperature phase transitions in the early Universe, with applications 
ranging from problems involving CP violation and baryogenesis, topological defect forma- 
tion, inflation, etc (for a short list of the different applications where SNR and ISB have 
been used see, e.g., Refs. |0-§[). 

ISB or SNR are direct consequences that in two or multi-field theories some of the 
coupling constants between fields can be negative while the model is still bounded from 
below. This can be the case in any extension of the standard model with a large scalar 
sector. Under these conditions, there can be an enhanced symmetry breaking effect at high 
temperatures, when thermal effects are taken into account in the model. This effect is clear 
when one considers a one-loop analysis of the effective potential in a simple model of two 
interacting scalar fields or, as in a 0{N) x 0{N) model. However, this simple analysis 
is too naive for two main reasons: first, at high temperatures the perturbative expansion 
in quantum field theory becomes, in most cases unreliable. This happens because there 
can be parameter regimes where powers of the coupling constants become surmounted by 
powers of the temperature, or due to the appearance of infrared divergences close to critical 
temperatures (as in field theories displaying a second order phase transition or a weakly 
first order transition). Second, the constrain conditions on the coupling constants, which 
may be important for the observation of ISB/SNR, usually requires large values for the 
couplings, in which case higher order loop corrections may become important as well and 
change appreciably the parameter space for ISB or SNR. 

In order to account for the above problems of the one-loop approximation, different 
methods have been used to analyze the question of ISB in quantum field theory at high 
temperatures. The results, however, have shown to be highly controversial, either by finding 
no evidence for ISB/SNR or favoring the phenomenon. Examples of the former were ob- 
tained in the context of methods like the large-N expansion [0, Gaussian effective potential 
chiral Lagrangian technique |^ and Monte Carlo simulations on the lattice |[T^. Some 



applications which find evidence for ISB/SNR are Refs. ||TT],|T2[, whose authors worked with 



self-consistent gap equations; Ref. |jT3|, also in the context of the large-N expansion, but 



reaching a different conclusion; Refs. [0,0 in the context of the renormalization group 



equations and Ref. |T^, also in the context of Monte Carlo simulations but this time sup- 
porting ISB/SNR. 

It is then obvious that it would be interesting to have a method to clarify this ques- 
tion without involving the possible difficulties related to the previous methods (like nu- 
merical precision in the Monte Carlo simulations, self-consistency in the gap equations or 
resummed perturbative methods, etc) used to study this problem. In this paper we use a 
nonperturbative technique known as the linear 5-expansion (also known as optimized per- 
turbation theory) (for earlier references see, e.g., [|l9l) and apply the method to the 
0{Ni) X 0{N2) scalar model to obtain the thermal masses to second order in the perturba- 
tive parameter 5. We then investigate their high temperature behavior to conclude about 
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the possibility of ISB/SNR. This model has already been considered before by Bimonte and 



Lozano in |]lT| where SNR was studied by means of the solutions of the gap equations of 
the model. By studying this theory we can also extend our results to the simple model 
consisting of two-interacting scalar fields, which has been extensively studied before. In 
this case the symmetry group corresponds to Z2 x Z2. This step will allow us to compare 
our results with the ones furnished by alternative methods. In a previous paper we 
have employed the optimized linear 5-expansion to study the resummation of higher and 
leading order thermal corrections showing that the use of a proper optimization scheme is 
equivalent to self-consistently solving the gap equation for the thermal mass, where lead- 
ing and higher order infrared regularizing contributions are nonperturbatively taken into 
account. An advantage of the linear ^-expansion is that the same simple propagator is used 
in the evaluation of any diagram, avoiding the potential bookkeeping problems associated 
to other resummation methods. Apart from being a powerful nonperturbative method, the 
optimized perturbation theory was originally formulated as a general theory applicable to 
arbitrary systems including strong interacting models, which makes it particularly interest- 
ing to use in connection with ISB/SNR where the issue of large coupling constants is an 
important one. 

This work is organized as follows. In Sec. II we introduce the model and the lowest order 
one- loop result. In the same section, the linear 5-expansion technique is briefly described. It 
is then used, in Sec. Ill, to evaluate the thermal masses up to order-5^ in the 3 + Id model of 
interacting scalar fields with global 0{N\) x 0{N2) symmetry. These calculations explicitly 
include two-loop momentum independent as well as momentum dependent diagrams with 
equal and different internal propagators. In Sec. IV we present our optimization results 
for the thermal masses and investigate, numerically, the possibility of ISB in the model. 
We specialize to the case A^^i = 90 and N2 = 24, where the model can be thought as 
representing the Kibble-Higgs sector of a SU{5) grand unified theory and also to the case 
A^^i = A^2 = Ij where it reduces to the Z2 x Z2 model, which has been the object of many 
studies in connection with ISB. Our predictions for the critical temperatures and size of 
the ISB parameter region are compared with the ones found in the literature. In Sec. V 
our concluding remarks are given. Two appendices are included to present some technical 
details and for a brief discussion of renormalization in the model. 



II. THE LINEAR (5-EXPANSION APPLIED TO THE EVALUATION OF THE 
THERMAL MASSES IN THE O(iVi) x 0{N2) MODEL 

A. The model and one-loop results 

In this work we consider the scalar 0{Ni) x 0(A^2) model described by 
2 



1=1 

where 



2 ^ 2i2 "^^ J,4 



A 



l<l>l + C,t, (2.1) 



1=1 



\c^\^l (2.2) 
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represents the counterterms needed to render the model finite. Note that Cct requires an 
extra piece if one attempts to evaluate the thermal effective potential which is not 
the case here. The boundness condition for the model described by ( |2.1| ) requires that the 
coupling constants satisfy the inequalities 



Ai > 0, A2 > and A1A2 > OA^ 



(2.3) 



As noticed by Weinberg this boundness condition allows for negative values of the 
cross coupling A which may lead to the nonrestoration of symmetries at high temperatures. 
Another possibility is that theories which are symmetric at T = may be broken at large T 
due to this negative value of A. The thermal masses for this model have been first calculated 
using the one-loop approximation which, at high T, gives 



Ml 



24 



and 



^ ~ TTin 



rp2 
24 



A.(^)+AiV. 



(2.4) 



(2.5) 



Let us consider the interesting case where A < and setting m\ and m| positive to have 
a symmetric theory at T = 0. Inverse symmetry breaking takes place if one chooses, for 
example 



A > 



No\ 3 



(2.6) 



which makes the coefficient of Ml negative while the same coefficient for M| is kept 
positive, due to the boundness condition. In this case, high temperatures will induce the 
breaking 0{Ni) x 0(A^2) 0{Ni — 1) x 0(A^2) at the critical temperature 



2 

mf 



24 



\N2 - Ai 



(2.7) 



The one loop approximation results will be investigated numerically and compared to our 
results in Section IV. 



B. The interpolated model 



The optimized linear (5-expansion is an alternative nonperturbative approximation which 
has been successfully used in a plethora of different problems in particle theory [ll^,pO|-|24| , 
quantum mechanics [p5| , p6| , statistical physics nuclear matter p8[ and lattice field 
theory One advantage of this method is that the selection and evaluation (including 
renormalization) of Feynman diagrams are done exactly as in ordinary perturbation theory 
using a very simple modified propagator which depends on an arbitrary mass parameter. 
Nonperturbative results are then obtained by fixing this parameter. The standard applica- 
tion of the linear 5-expansion to a theory described by a Lagrangian density C starts with 
an interpolation defined by 
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C' = {1- 5)C^{'n) + 5C = Coiv) + S[C - A 



(2.8) 



where CqIt]) is the Lagrangian density of a solvable theory which can contain arbitrary mass 
parameters (rj). The Lagrangian density interpolates between the solvable Cq{vi) (when 
5 = 0) and the original C (when 5 = 1). For the present model one may choose 



^ 2/2 



2^. 



2 J,2 



(2.9) 



and following the general prescription one can write 

5 



=1 L 



4^ 



4 



(2.10) 



where Vtf 
model finite, is given by 



The term >C^,., which contains the counterterms needed to render the 



E 



4! 



5C 



+ -5i?f(r^i,r^2)0.' 
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6C'<p\<Pl, (2.11) 



where Af, 5f, Cf and are the counterterms coefficients. One should note that the 
5-expansion interpolation introduces only "new" quadratic terms not altering the renormal- 
izability of the original theory. That is, the counterterms contained in as well as in the 
original Let-, have the same polynomial structure. 

The general way the method works becomes clear by looking at the Feynman rules gener- 
ated by . First, the original 0f vertex has its original Feynman rule modified to —iSXi 
(the same applies to the mixed vertex (flcf^)- This minor modification is just a reminder 
that one is really expanding in orders of the artificial parameter 6. Most importantly, let 
us look at the modifications implied by the addition of the arbitrary quadratic part. The 
original bare propagator. 



S{k) = i{e - + ley 



becomes 



s{k) = i{e - + ze)-^ 



k"^ — + ie 



1 - 



k^ 



mf + ie 



(2.12) 



(2.13) 



indicating that the term proportional to T]f(j)f contained in Co is entering the theory in a 
nonperturbative way. On the other hand, the piece proportional to Srjfcfi is only being 
treated perturbatively as a quadratic vertex (of weight iSrjf). Since only an infinite order 
calculation would be able to compensate for the infinite number of {—irjf) insertions con- 
tained in Eq. (|2.13|) , one always ends up with a rji dependence in any quantity calculated 
to finite order in 6. Then, at the end of the calculation one sets 6 = 1 (the value at which 
the original theory is retrieved) and fixes rji with the variational procedure known as the 
Principle of Minimal Sensitivity (PMS) fST 







(2.14) 
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where P represents a physical quantity calculated perturbatively in powers of S and then 
extremized with relation to the r]i parameters^. This optimization procedure applied to the 
thermal masses will be discussed in Sec. IV. 



III. THE THERMAL MASSES UP TO ORDER-J^ 

We can now start our evaluation of the thermal masses, defined by 



(3.1) 



where (p) is the thermal self-energy. At lowest order (first order in S) the relevant contri- 
butions, which are momentum independent, are given by = 1,2 and i ^ j) 



2 V 3 J Jt {27tY k'^ -VLi + it 2 ^ Jt {2'nY k'^-VL] + it 



d'^k 



■ (3.2) 



The temperature dependence can be readily obtained by using the standard imaginary time 
formalism prescription 



Then, the self-energy becomes 
Si (p) = ~Sr]^ + ST— 



It {2nY^''^Vl 



{2ttY-^ 



(3.3) 



E 



d'^'^k 



{2iTr~^ ul + E, 



(2vr)^-i ul + Ey 
(3.4) 



where E"^ = \s? + VL^. Summing over Matsubara's frequencies one gets 

X fNi + 2\ f d^-^k f 1 1 



(p) = -Sv^ + 6 
.A 



2V 3 J J {27tY-^ \2Ei Ei[l - exp{E/T)] 
d'^-^k f 1 1 1 



{2nY-^ \2E, E,[l-eMEj/T)]j ' 
Then, using dimensional regularization |^ (rf = 4 — 2e) one obtains the thermal mass 



(3.5) 



327r2 V 3 



A 



+ 327T^T^h i ^ 



+ 32'K^T'^h ( ^ 



(3.6) 



^For a discussion of convergence as well as renormalization in the method, see Ref. and 
references therein. 



5 



where /x is a mass scale introduced by dimensional regularization and 



oo ^2 



47r^ Jo [x2 + 2 [exp(x2 + yf)2 - 1] 

Note that the temperature independent term diverges and must be renormalized. In this 
paper we chose the Minimal Subtraction (MS) scheme where the counterterms eliminate the 
poles only. At this order the only divergence in ( p.6|) is 



which is easily eliminated by the 0(5) mass counterterm 

= Bt\n,,n,) = J— + Aiv,fi^ ) . (3.9) 



By looking at Eq. (|3.6| ) one can see that the terms proportional to 5 represent exactly those 
that appear at first order in the coupling constants in ordinary perturbation theory excepted 
that we now have Vtf instead of mf, 5Aj instead of Aj and 5A instead of A. Therefore, it is 
not surprising that to this order the renormalization procedure implied by the interpolated 
theory is identical to the procedure implied by the original theory at first order in the 
coupling constants. 

Let us now analyze the temperature dependent integral which is expressed, in the high 
temperature limit (?/j = f2j/T ^ 1), as 



him) = YA-hy'- T^^^' 



In ( ^\ +1E-]: 



+ 0{y^) . (3.10) 



Art:, 

In principle, since rji is arbitrary, one could be reluctant in taking the limit Vti/T <ti 1. 
However, as discussed in our previous work |^0[, the use of both forms for the integral h{yi) 
does not lead to any significant numerical changes in the optimization procedure. This is 
also true in the present work. Then, by taking h{yi) in the high temperature limit, one 
obtains the 0{6) thermal mass: 

M^ = n^- 5r,^ + 5\ (^^) X,{T) + 5\N,X,{T) + 0{5') , (3.11) 

where we have defined the quantity Xi{T) 

X,(T) = — - + -^L(T) , (3.12) 



with L{T) given by 



L(T) = In f—^j - 7i. . (3.13) 



For notational convenience when expressing the remaining contributions to the self-energies, 
we also define the additional quantities Yi(T), Zi{0), Wi{0) and Ri(T) given by 
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(3.14) 



In 



47r/x2 



+ 7e 



2 2 

+ 12' 



(3.15) 



iy,(o) 



In 



47r/i2 



1 

+ 2 + 12' 



(3.16) 



and 



(3.17) 



At (9(5^) the sclf-cncrgy receives contributions from momentum independent as well as 
momentum dependent diagrams. Let us first consider the momentum independent diagrams, 
which to this order have one and two loops. The order-(5^, one loop, momentum independent 
contribution is given (in the high temperature limit) by 



r2. 2iNi + 2 



1 



327r2e 



(3.18) 



and 



-52 

-'i,3 



327r2e 



(3.19) 



As discussed in Appendix B, the above contributions can be rendered finite using the mass 
type counterterms contained in C^^ which are tailored to account for divergences arising from 
the extra quadratic vertices introduced during the interpolation process. The momentum 
independent two loop contribution is given by the four "double scoop" diagrams 



/TV- + 2 



+ 



m 

(167r)2 



1 



(327r2)2e2 



8X,{T) 



27r 



+ 



327r2e 



[X,{T)+nT)] 



(327r2)^ 



[Zi(0) + Wim 



J3 



3847rf],- 1287r2 



(3.20) 



{Ni + 2) 



^1 



(327r2)2e2 327r^€ 



[Xj{T) + R,{T)] 



+ 



m 

(167r)2 



+ 



(327r2)^ 



[Z,{Q) + WM]} , 



+ 



3847rQ,- 1287r2 Q,- 



(3.21) 



{N, + 2) 



[X,{T) + Y,{T)] 



+ 



L{T) 
(167r)2 



. ' 27r 



(327r2)2e2 327r2e 



(327r2)2 



rp3 



+ 



3847rfij 1287r2 



(3.22) 
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and 



[X,{T) + R,{T)] 



SSAnilj 1287r2 ilj 



+ 



L{T) 
(167r)2 



8X,(T) - 



+ 



27rfijJ (327r 



2^2 



[^.(o) + w^.(o)] . 



(3.23) 



To render these diagrams finite one needs mass and vertex counterterms |31[]. Considering 



the 0{5) mass counterterms used to ehminate the divergences in S^lvi (see Eq. 
is able to build two one loop 0{5'^) diagrams whose contributions are given by 



_ 5^\{N, + 2) 



and 



■'1,9 



3(327r2)2 



2^2 



(327r 



A,-fi|(Ar, + 2) 



+ An,' AT, 



32_7r2 



)), one 
(3.24) 

(3.25) 



Additionally, from the vertex counterterms appearing in ( ^.11| ), with Cj and C (at order 
(9(5^)) given by 



5^ 



327r2e 



and 



o^2\2 f2\ 

C = + [Ai(Ari + 2) + A2(Ar2 + 2)] 



one can built the one- loop (9(5^) diagrams: 



.52 



(A^^ + 2) 
167r2 3 



A 



and 



-'Lll 



6 



327r2 
327r2e2 



+ A 



2 3 AT,. 



1 Q2 

327r2e2 e ^ ^ 327r2 ' 



(3.26) 



(3.27) 



(3.28) 



4A2 + ^ [A, (AT, + 2) + A, (AT, + 2)] 



- -X,{T) + 



327r: 



rWG(0) 



(3.29) 



The next contribution to the self-energy at 0(5') comes from the two- loop "setting sun" 
diagrams. These momentum dependent contributions are given by the setting sun diagram 
with equal mass internal propagators, 



-6 



+ Giiifi) 



(3.30) 



and by the one with internal propagators with different masses. 
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(3.31) 



In the above expressions, Gq is the zero temperature part of the diagrams and Gi and G2 
are the finite temperature ones (with one and two Bose factors, respectively). The details 
of the evaluation of ( p.30| ) and ( |3.31|) are given in Appendix A. Their contributions are 



3(327r2)2 

^2/ 



p2 WK^Xi{T) 



6(47r) 



+ 5 



2\'(A^. + 2) 



In 



1^? 



+ 5' 



3(47r)2 \^47r/i2 

24 ~ "8^ ~ 16^ 
A2(Ar, + 2)T2 " 



In 



UttT/ 



2 + IE 
+ lE 



72(47r)^ 



In j + 5.0669 



1^? 



471/^2 



1.9785 



(3.32) 



and 



2 (47r) 



n 



In(r2^) ) - 



+6' 
-52 
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(47r)4 
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(47r)4 


A^AT, 





1 + 



2 (47r)4 
^2^2 r 
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A^l 

(47r)2 7 



1 / 



3 - 27ij - 2 In 



7 + ^ - 67i5 + 27I - 2(3 - 27s) In 
6 



n 



1 H - H ln(?^ 



1 

-■i 
2 



-n^ In^(n^) + 



'1 



21n^ 



47r/i2 



TT 



Li2(l-n^)- — 



3 - 27s - 2 In 



47r;u2 



(47r) 



In 



^47r/i2 



2 + 7i? 



+ ln(ra^ 



In 



^47r/i2 



2 + 7i? 



2m2 r 



+5' AT, 



2 In 



^^2, I 2^^j^ 

ST 



+ 5.0669 



(3.33) 



where n = Qi/Qj. 

From Eqs. (pj), (|3l8|)-(|3:25|), (pD, dOl , dpl) and (lOl one easily sees that all the 
temperature dependent divergences cancel exactly. The remaining divergences are handled in 
the usual way, being canceled by the counterterms appearing in Cct, Eq. (|2.2| ) (see Appendix 
B). The sum of the remaining finite terms of each contribution to the self-energy makes the 
total contribution to the thermal mass up to order 5^, Mf = Qf — 5r]f + Sj^i + . . . + Sj_i3. 
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IV. NUMERICAL RESULTS 



We are now in position to set 5 = 1 and apply the PMS to the thermal masses. Before 
doing that few points concerning the optimization procedure should be clarified. First, 
we recall that our interpolation procedure has been carried out in a very general way by 
assigning a different interpolation parameter to each field. Although general, this procedure 
brings in two arbitrary parameters which have to be fixed with the PMS. Of course this 
complicates the numerical optimization procedure, which is needed at higher orders, since 
one has to look for extrema in the r^i, rj2 space. Sometimes these extrema show up as saddle 
points which are hard to detect numerically. In principle, bearing in mind that 0i and 02 
are both the same type of fields, one could be tempted to use the freedom allowed by the 
interpolating process to set rji = ri2 = r]. In order to assess the validity of such choice we 
have done the optimization in both ways finding that the case Tji ^ ■r]2 gives better results. 
The second point regards the actual quantity to be extremized. By looking at our equations 
one can easily see that, except for the case where Ni = N2, ml = m\ and Ai = A2, the 
PMS applied separately to Ml and M|, at the same temperature, can generate different 
values for the same f]i. One way to avoid this would be to evaluate and optimize a more 
comprehensive quantity such as the effective potential. In fact this procedure has been 
advocated in other applications of the optimized 5-expansion However, in the present 
work we have not attempted to evaluate the effective potential since, at 0{5^), this would 
imply in the evaluation of three-loop zero point functions which become rather cumbersome 
at finite temperatures 0. Although this may be regarded as a controversial point within 
the method it will not be discussed any further in the present application. Here we follow 



the original prescription given in where it is suggested that the PMS should be applied 
to each different physical quantity so that rji can be adjusted to the relevant energy scale. 
The validity of such procedure will be judged by comparing our results with well known 
predictions in the A^i = A^2 = 1 limit. Finally, as discussed in [^], one can easily see that 
applying the PMS to at first order in 5 produces coupling independent values for fji 
which do not furnish truly nonperturbative results. Therefore we will only investigate the 
results generated at order-5^. 

We work in units of the arbitrary scale /i introduced by dimensional regularization, by 
extremizing the dimensionless quantity Mf/ix^. To search for inverse symmetry breaking 
at high temperatures one must ensure that the symmetries are restored at T = 0. This is 
achieved by setting mf > and by observing the boundness condition while setting A — *• —A 
in all our equations. 



A. The Z2 X Z2 case. 

Let us start by studying the A^i = ^"2 = 1 case where the theory reduces to the Z2 x Z2 



model which has been extensively studied in the literature. In (see also [|T5|) it is 



claimed that the critical temperature is not reliably estimated by the loop expansion which 



See Ref. [33| for a discussion concerning the J-expansion evaluation of the effective potential. 
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does not include temperature effects in the coupling constants. On the other hand, the 
nonperturbative renormalization group approach (RGA) used in 0] does include those 
effects and as a first check of the validity of our results we estimate the critical temperature 
for the Z2 X model comparing our results at (9(5^) with the ones furnished by the two 
methods mentioned above. This is done in the first two rows of Table I and also in Fig. 
1. In all cases m\j ^ = ml/fi"^ = 1.0. Just above the critical temperatures presented in 
the first two rows the thermal mass is negative whereas M| is positive. Our numerical 
optimization strategy is the following: for a temperature below the value predicted by the 
1-loop approximation (ILA) we identify the extremum in the rji, 772 space. Both masses in 
the extremum positions are found to be consistently positive and the extremum positions 
are found to be unique (and consistent with the high temperature approximation, in the 
sense that < T). The extremum is then followed as the temperature is increased. When 
one of the masses becomes negative the value of is obtained and the values of f/i, f/2 are 
registered. For example, around the critical temperature the values obtained by extremizing 
Mf are f/i ^ 0, fj2 ^ 0.25 T and f/i ^ 0, f72 ~ 0.15 T for the parameters in the first and 
second rows, respectively, in Table I. From this table one can see the remarkable agreement 
between our results and the ones provided by the RGA. In Fig. 1 we compare the numerical 
values of the critical temperature with the ones predicted by the ILA and RGA for several 
values of the cross coupling. Once again our results agree with the ones given by the RGA 



Tj,|T5|. We have also investigated the behavior of the critical temperature as a function of A2 



for fixed ml 2/ /^^ = 1-0, Ai = 0.009 and A = —0.025. As can be seen directly from Eq. ( p. 41) , 
the one-loop predicts that Tc/rrii ~ 49.0 for any value of A2, while the RGA predicts that 
it increases with the latter. Once again our results, not shown, agree with the ones given 
by the RGA. These findings give us confidence about the correctness of the optimization 
procedure adopted. 

The critical temperatures predicted by the 5-expansion, as well as by the RGA, are 
higher than those predicted by the ILA. This could be roughly understood by recalling that 
for the phase transitions considered the PMS applied to the mass which signals the phase 
transition (M^ for the parameters of Table I) generates a f/2 which is an increasing function 
of T whereas it gives f/i ^ 0. Therefore, at high T, the optimized mass 02 which "dresses" 
the 5-expansion propagator associated with 02 is greater than its counterpart Clf. In M^, 
the quantity fl"^ is always associated with the coupling A which drives ISB. It seems that the 
increasing of f^l with the temperature suppresses the effect of A in enhancing ISB. Another 
interesting point to be discussed regards the inclusion of temperature dependent coupling 
constants. While this possible dependence is overlooked in the ILA it is taken into account 
in the RGA. These effects are more subtle to be observed directly in our case. However, 
by looking at the nature of the diagrams considered by us at order-(5^, one may have an 
idea of how these effects enter our calculations. Basically, one can think of double scoops 
and setting suns diagrams as being tadpole diagrams with vertex corrections. Expanding 
the vertex to order-5 gives the simple one loop diagrams considered in Eq. ( |3.11|) , while 
expanding the vertex to order-5^ gives the double scoops and setting suns which indirectly 
contain one-loop temperature dependent corrections to the vertex. 
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B. The 0{Ni) x 0{N2) case. 



Let us now consider the more realistic 0{Ni) x 0{N2) case. We follow Bimonte and 
Lozano, in the first reference in by choosing A^^i = 90 and N2 = 24 so that the model 



can be thought of as representing the Kibble-Higgs sector of a SU{5) grand unified model. 
Let us start by estimating some critical temperatures and comparing our results with the 
lowest order one-loop results given by Eqs. (U) and (U). The two first rows of Table 
II display our results for the set of parameters considered in the previous subsection. Our 
predicted values for the critical temperatures are proportionally much higher than for the 
A^^i = A^2 = 1 case. The two last rows of Table II display the results for more typical 
parameter values for this model, which are usually taken around Ai = 0.8 and A2 = 1.0 (see 
1^). These results suggest that the parameter region for ISB predicted by the 5-expansion 
is smaller than the one predicted by the lowest order ILA prediction. To illustrate that 
we offer Figs. 2 and 3. Figure 2 shows our (9(5^) result (dot-dashed line) compared with 
the results produced by the ILA at lowest order (thin continuous line) for fixed Ai = 0.8 
and T = 5.0 fi. The upper parabola represents the limiting region for boundness. The 
dashed line is the Bimonte and Lozano's next to leading order one-loop (BLA) result for an 
arbitrarily large T. Fig. 3 is similar to Fig. 2 except that Ai varies while A2 is kept fixed at 
the unity value. Note that the next to leading order correction considered in |TT] for the one- 
loop approximation reduces significantly the region of ISB given at lowest order. However, 
the next to leading order calculation of the BLA approach considers up to the second term, 
which is mass dependent, in the expansion of the temperature dependent integral given by 
Eq. (|3.7|) . This means that their procedure does not take into account logarithmic terms in 
the expansion of h{yi) since they would give corrections of order Xj In Aj which also arise from 
two-loop diagrams. Our calculations, on the other hand, avoid these bookkeeping problems 
and all contributions are consistently and explicitly considered in a given order. Moreover, 
support for ISB/SNR also arises in the work of Amelino-Camelia Ref. who, like us, 
considers two-loop diagrams in the simple Z2 x Z2 model. 

To illustrate how the temperature affects the ISB region we offer Fig. 4 where we have 
fixed Ai = 0.8 while A2 varies. It is clear that the ISB parameter region for T = 5.0 fi is 
smaller than for T = 10.0 fi. That is, increasing temperatures favor ISB. The parameter 
regions predicted by the ILA have the same temperature dependence although the difference 
in size is less significant than the one displayed in Fig. 3 for the 5-expansion. 

Finally, let us investigate the possible patterns for ISB. According to the one- loop ap- 
proach (to lowest and next to leading order) as well as to the RGA, there are only two 
possible phases at high T: either the theory is completely symmetric or one of the two 
symmetries is broken. However, including two-loop contributions and going beyond the 
simple perturbative expansion could alter this picture introducing a third possibility where 
the two symmetries are broken at high T. To analyze this possibility one must increase the 
temperature beyond the critical values shown in Tables I and II which display only the 
connected with the phase transition from the symmetric phase to one of the nonsymmetric 
phases. Surprisingly we find, that for some parameter values, the mass associated with the 
symmetry which survives the first transition has a tendency to decrease as the temperature 
increases beyond the first critical value. As an example of this, we show in Fig. 5 the results 
for the case analyzed in the first row of Table II for both and M|. M| is positive right 
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above the critical temperature shown (now labeled Td) which is associated with the break- 
ing in the 0i direction. We sec that, very quickly, M| becomes monotonically decreasing 
with the temperature and becomes negative, through a second order phase transition, at 
Tc2/^2 — 18.9. The same behavior is observed for the parameters shown in the second 
row of Table II, where we observe symmetry breaking in 02 direction at Tc2/m2 — 139. 
However, analyzing the other two cases shown in the third and fourth rows of Table II, we 
find that the would be mass associated with the second symmetry breaking increases mono- 
tonically with the temperature, signalling that this symmetry possibly remains unbroken 
at high temperatures. A similar behavior happens for the two cases analyzed in Table I, 
for Ni — N2 — I. Although M| initially shows a decrease in value as the temperature is 
raised, it soon becomes monotonically increasing with T. Note that just by interchanging 
the values of Ai and A2 one trivially observes yet another pattern this time with Tc2 < Td. 
Which pattern will be actually followed depends on our initial choice for the values of the 
couplings. 

Therefore, our results suggest that a possible symmetry breaking along the second field 
direction takes place, for large values of N, in a narrow region of parameters. It is possible 
that this alternative symmetry breaking pattern will show up only in nonperturbative cal- 
culations which consider up to two loop terms. It would be interesting to investigate this 
possibility using other nonperturbative approaches. 

V. CONCLUSIONS 

We have used the optimized linear S expansion to investigate inverse symmetry breaking 
at high temperatures using multi-field theories. Our order-5^ calculations take full consid- 
eration of two-loop contributions, including the momentum dependent "setting sun" type 
of diagrams. To our knowledge, a complete calculation associated with the phenomenon 
of ISB/SNR which includes these contributions in the 0{Ni) x 0(A^2) model has not been 
fully considered before. In order to asure the reliability of the method we have started with 
the scalar Z2 x Z2 model which has been extensively investigated in connection with inverse 
symmetry breaking problem. We have shown that our optimized results agree well with 
those obtained with the Renormalization Group approach, especially as far as the critical 
temperatures are concerned. This has allowed us to establish the 5-expansion as a reliable 
nonperturbative technique to investigate ISB. We have then investigated the more realistic 
scalar 0{Ni) x 0(A^2) model which may be related to the Kibble-Higgs sector of a SU{5) 
grand unified model. All our results strongly support the possibility of inverse symmetry 
breaking (or symmetry nonrestoration) at high temperatures. Surprisingly, we have also 
found evidence for a second phase transition taking place for some values of the couplings 
and large values of N. According to our results two other possible high temperature inverse 
symmetry breaking patterns are 

0{Ni) X 0(N2) ^ 0(Ni - 1) X 0(^-2) ^ 0(Ni - 1) x 0(N2 - 1) , 
where Td < Tc2, or 

O(iVi) X 0{N2) ^ 0{Ni) X 0{N2 - 1) ^ 0{Ni - 1) x 0{N2 - 1) , 
where > Tc2- 
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APPENDIX A: 



Consider a general setting sun diagram given by = 4 — 2e): 



dfiq 



■ (Al) 



At finite temperature we express the momentum integrals as in ( |3.3| ). The discrete sums in 
the Matsubara frequency can easily be done if we reexpress Eq. ( |Al| ) in terms of the Fourier 
transformed, in Euclidean time, expressions for the field propagators. 



^ , . , _ 4e /• d^-^kl d^-'k2 d^-^h 
53(p-ki-k2-k3) , 



(A2) 



where Gj(k, r) is the propagator, which can be written in terms of a spectral function p(k, r/), 
as 



G.(k,r)= r°'^[l+n(r/)]p,(k,r/)e-^l^l 

J-oo ZTT 



where n is the Bose distribution and 



(A3) 



(A4) 



Using ([A3D in (|A^) , the identities e*^"^" = 1, n{ri) = e"^''[l + ^(77)], 72(77) = -[l + '^l-^)] and 
performing the r integration, we get 



Gijj(p,iiu;„)=/i'^y" 



(2 



TT 



■27r 



(2 



TT 



2n 2n 2n 



■Pi(ki, ?7i)pj(k2, r]2)pj{ks, 773) 



X [l + n(r7i)] [l + n(r/2)] [1 + ^(r/g)] 
X(53(p-ki -ka - kg) . 



1 



1 



ii^n + m + ^2 + r]3 ivn - ^ - r]2 - r]^ 



(A5) 



Performing the analytic continuation in the above expression, iun —>■ Po + i^, and using 

e 



lim 



7r5(a;) 



(A6) 



ReGijj + ilmGijj. The 



we can separate Gijj in real and imaginary contributions, Gijj 
imaginary contribution gives the field decay rate and it is not important in the present 
work. The real part contributes to the thermal mass and it is given by (after performing 
the /c3 momentum integral with the help of the Dirac delta-function in (|A5|)) 
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ReGijj{p,po) = J 



4e 



{27rY~^ (27r)'^-i 8E,{k,)E,{k2)E,{p - ki - k^) 
X [6ivi - E,ik,)) - 5(7], + E,(ki))] [6iv2 - E^ik^)) - 5ir]2 + E^ik^))] 
X [5(r73 - E,ip - ki - k2)) - (5(r73 + ^.(p - ki - k^))] 

^ ^ ' - ki - ka - kg 



dr]idr]sdr]-i 



(A7) 

_P0 + Vl+V2+V3 P0-V1-V2-V3. 

After performing the r] integrals with the help of the Dirac delta-functions and after some 



algebra, Eq. (|A7|) can be written in terms of three terms: a temperature independent one, 
which is just the zero temperature contribution Gyj^, and two other terms with one and 
two Bose factors, which gives the Gjjj,! and Gjjj,2 terms appearing in Eq. (|3.31| ). 

The zero temperature contributions in both cases have already been evaluated in the 



literature. The contribution Gnj,o(p) has been evaluated in details in p4| , p5| where the 
quoted result for the on mass shell {p^ = Qf) case is 



Giiifl{p) 

which gives 

G Hi flip) = - 



(47r2)4(l-e)(l-2e) 



2(47r)4 

3a' 



'47r/i^ 



3 1 
"2? ^4^ 



3 - 27£; 2 
m 



+ 



19 



4(47r)4e 



(4vr) 



In^ 



1^? 



47r/i2 



+ 1.9785 



(A8) 



(A9) 



where we purposefully left the momentum dependence in the relevant divergent term to 
make explicit the need for a wave-function renormalization counterterm. 

The zero temperature contribution to the mixed setting sun diagram is also given in Ref. 
for the on shell case (p^ = Qj) can be written as 



Gijjflip) 



+ 



11 



(4vr)Ml-e)(l-2e) 



n 



n 



2e2 



2 rl 

- + ln(n2) 



1 + -^n' -(1 + ^1 \n{n^) - -n^ \n^{7?) + 



'1 — 77, 



2\2 



Li.fl 



n — 



71 

~6 



(AlO) 



where n ■ 

Gijjflip) 



Vli/Vlj and 112(2:) = YaLi /l'^- We then obtain the result 



(47r)4 [e2 V 

4?^-T^"^"' 

(4vr)4 1^+2] 



n 

y 



1 / 



3 - 27b - 2 In 



7+ — ~Q^e + 2Ye- 2(3 - 2^e) In 
o 



+ 



(4vr)' 



11 / tt^x 1 

2 / 1 I 1 1 / 2\ 2 1 2/ 2n 



47r/i2 



2 In' 



47r;u2 
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2\2 



(1-^2) 



Li2(l - n") - 



TT 



(4vr) 



3-27^ 



21n ' ^ 



(All) 



The finite temperature terms Gjjj.i and 6*^^.2 can be worked out as follows. Taking p = 
in these terms, allows us to reexpress them by (using ki.k2 = kik2 cosO) 



ReG',,,-i(0,po) 

d 



d'^-^ki d'^-^k2 



1 



X 



n(E,(ki)) 



(27r)^-i (27r)'^-i 8Ei{ki)Ej{k2) kik2 
In { [pI - (E,(ki) + E,(k2) + E,ik^ + k2)f 

O cos (7 ^ L 



-E,(ki) + E,{k2) + ^,(ki + k2)f] } 



+2n(E,(k2))— — In { [pI - (E,(ki) + E,(k2) + i?,(ki + k2)f 



p2 - (E,(ki) - i?,(k2) + i?,(ki + k2)f] } 



(A12) 



and 



ReGjj,,-2(0,po) = At 



4e 



1 



X 



(27r)^-i (27r)'^-i 8Ei(ki)E, (k2) A;ifc2 
X [2n(E,(ki))n(E,(k2)) + n(E,-(ki))n(E,(ki + ks))] 

In { [pI - (E,(kO + E,(k,) + E,(ki + k,))^ 

-E,(ki) + E,(k2) + ^,(ki + k2))'" 
pI - (E,(ki) - i?,(k2) + i?,(ki + k^))'^ 
pI - (Ei(ki) + E,{k2) - ^,(ki + k2))'] } 



X [vl 



(A13) 



For i = j (equal mass propagators) Eqs. ( |A12|) and (|A13|) give the same expressions ob- 
tained by Parwani in Ref. In special, we note that, in Eq. (|A12| ), the terms given by 
9/9cos^^ln{. . .} 2ki/k2 as ^2 oo. We can subtract and add this term in the appropriate 
places in ( |A12| ), obtaining the analogous expressions given by Parwani, 



ReGijj_i(0,po) — Fijjfi + Fijj^i + Fijj^2{Po) 



where 



and 



ijj,0 



(4^7 



(A14) 



(A15) 



E. 



(4vr) 



In 



T 



47r/i2 



2 -7s 



2 -7s 



(A16) 
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where in the above equations, h{yi) is given by ( |3.10|) . The remaining terms, -Fijj 2 and G'ijj,2 
can be evaluated on-shell and a similar contribution have already been computed earlier in 



the literature, see Ref . , from where we obtain 



2rp2 r 



3(4^ 



2 In 



5.0669 



(A17) 



The finite temperature contributions Gm^i and GjM,2 are given, as in Ref. and, from 
the previous equations, Eqs. (|A14|) - ( |A17|) , they can be written as 



6- 



Re[Gm,i(0, rij)] — Fiiifl + Fiii^i + Fi 



iii,2 



(A18) 



where0 



Finn — —5 



2A2(Ar, + 2)T2l^ /a 



3(47r) 



T 



(A19) 



and 



\1{N, + 2)T\ fVk 
T 



FiiiA — —5"^ .J — h 



3(47r) 



In 



47r/i2 



Fill ^2 (^j 



-Re[G'2(0,a)] 



72(47r) 



+ 2-7E 



In (^] + 5.0669 



(A20) 



(A21) 



Putting all these contributions together, we obtain the results given in Eqs. ( p3.32D and 
(PI). 



APPENDIX B: 



To obtain the total finite order 5^ contribution one can add all divergences appearing in 
Eqs. (gl8D-(^:25l), (|]29D, and (pl). As it can be easily seem all the temper- 

ature dependent divergences cancel exactly and one is left with temperature independent 
poles. 

By looking at all terms which contribute to this order one can identify two classes. The 
first is composed by diagrams such as the ones described by Eqs. (|3 . 20|) , (|3l21| ) , ( P . 22]) , ( P . 23| ) , 
(F^) , (^:25D , (lOSi ) , (|3:29| ) , i^M) and (pl). All of them are analogous to the diagrams 
which appear at second order in the couplings in the original theory and can be rendered 
finite by similar mass and wave-function counterterms. This procedure has already been 
illustrated at order-5. One can generalize this procedure by stating that diagrams belonging 
to a general order-^"-, and containing any combination A™'A^A2 such that m + k + I = n 



^Note that there is a misprint in Eq. (2.32) of Ref. |2^, where there is an extra 1/2 multiplying 
that equation. 
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will be renormalized exactly as when ordinary perturbation theory is applied to the original 
theory. Then, one just has to replace the original masses with the relevant interpolating 
masses f^i 2. It is easy to check that for those diagrams the most divergent terms will display 
e~" poles. 

The second kind of diagram is exclusive of the interpolated theory and carries at least one 
5r(f (or 5r(^) vertex. At (9(5^) these diagrams are described by Eqs. ( |3.18|) and ( p.l9| ), which 
display the divergent rjf and rj^ pieces. Looking at one identifies a r/j^jT^J-dependent 
coefficient whose Feynman rule is (?7i, 772)- Since the actual pole is of order-5^ one 
identifies this coefficient as having the same structure as the mass counterterm (fii, f22), 
displayed in Eq. (|3.6|), except that we now have rlf - instead of Vtf -. Therefore, (9(5") 
diagrams belonging to the second class will make use of the counterterm 55"^" (771, 772). This 
coefficient is similar to "(f2i,f22) which has been evaluated in a previous order. One 
can also easily check that for these diagrams the most divergent terms will have e"~^ poles. 
Moreover, power counting reveals that those 5rjf (or 5rij) insertions make the loops more 
convergent. For example, all one loop diagrams of order 0(5"), with n > 3 are finite. 



The renormalization prescription adopted here is analogous to the one shown in for 
the one-field case, where we have shown that the order by order renormalization holds at 
any higher orders in 5. In Ref. EO] the renormalization procedure is treated in more detail 
for the simple case. 
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TABLES 

TABLE I. Results for the critical temperature as obtained in the 1-loop approximation, the 
renormalization group approach in Ref. [14] and in the (5-expansion, respectively. 



Ai 


A2 


A 








0.06 


1.8 


- 0.1 


24.5 


33.6 


33.3 


0.0167 


0.6 


- 0.025 


53.7 


67.9 


66.9 



TABLE II. Results for the critical temperature in the 0(90) x 0(24) scalar model. 



Ai 


A2 


A 




rf /mi 


0.06 


1.8 


- 0.1 


6.5 


15.3 


0.0167 


0.6 


- 0.025 


16.5 


65.4 


0.9 


1.0 


- 0.141 


2.2 


5.0 


0.8 


0.7 


- 0.091 


3.4 


10.0 
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Figure Captions 



Figure 1: The critical temperature Tc/mi as a function of A (< 0) for A^^i = 1 and 
= 1 with the following parameter values: m^g//^^ = l-O^ -^i = 0.018 and A2 = 0.6. The 

dashed line is the one-loop prediction, the continuous hue is our result and the dots represent 

the values obtained with the RGA. 



Figure 2: Region of ISB at fixed Ai = 0.8 for Ni = 90 and N2 = 24. The thick parabola 
limits the region for which the potential is bounded. The continuous and dashed lines are 
the zeroth order (at T/fi = 5.0) and first order (at an arbitrarily large T) results given in 



Ref. |Tl|. The dot-dashed line is our result at T/fi = 5.0. The region of ISB is the one in 



between the boundness curve and the other curves for each case. 



Figure 3: Same as in Fig. 2, but now at fixed A2 = 1.0. 



Figure 4: The region of ISB for two different temperatures (Ai = 0.8): T/fi = 5.0 
(dashed line) and T/fi = 10.0 (thin full line). 

Figure 5: The behavior of (dashed line) and M| (upper full line) as a function of 
the temperature for Ai = 0.06 A2 = 1.8 and A = -0.1 (A^i = 90 and N2 = 24). 
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